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Abstract. Cells in multicellular organisms adhere to the extracellular matrix through two-dimensional 
clusters spanning a size range from very few to thousands of adhesion bonds. For many common receptor- 
ligand systems, the ligands are tethered to a surface via polymeric spacers with finite binding range, thus 
adhesion cluster stability crucially depends on receptor-ligand distance. We introduce a one-step master 
equation which incorporates the effect of cooperative binding through a finite number of polymeric ligand 
tethers. We also derive Fokker-Planck and mean field equations as continuum limits of the master equation. 
Polymers are modeled either as harmonic springs or as worm-like chains. In both cases, we find bistability 
between bound and unbound states for intermediate values of receptor-ligand distance and calculate the 
corresponding switching times. For small cluster sizes, stochastic effects destabilize the clusters at large 
separation, as shown by a detailed analysis of the stochastic potential resulting from the Fokker-Planck 
equation. 

PACS. 87.15.-v Biomolecules: structure and physical properties - 05.10.Gg Stochastic analysis methods 



1 Introduction 

Cells adhere to the extracellular matrix or to each other 
through a multitude of weak interactions. The transient 
character of their adhesions allows cells to adapt quickly 
to changes in their environment. In particular, cell adhe- 
sion has to be transiently down-regulated during impor- 
tant physiological processes like migration or division. For 
tissue cells, the main receptors for cell-matrix adhesion 
are integrins, which are linked on the extracellular side to 
ligands like fibronectin and on the cytoplasmic side to the 
actin cytoskeleton [1 . This provides structural integrity 
between extracellular matrix and cytoskeleton, which is 
important because cell-matrix adhesion clusters usually 
have to operate under considerable mechanical load. The 
behaviour of single receptor-ligand bonds under force was 
first discussed by Bell [2], who suggested that bond life- 
time is reduced exponentially by an applied load. This 
concept has been impressively confirmed and extended 
by dynamic force spectroscopy [3UH[5j[6l[7] , which showed 
that binding strength can only be defined in a dynamical 
context. In particular, it has been shown both theoreti- 
cally and experimentally that for linearly increasing load, 
binding strength increases linearly with the logarithm of 
loading rate. If loading occurs through a soft polymeric 
linker, binding strength is decreased compared to load- 
ing at the same speed but through a rigid linker [5]. In 
practise, cell adhesion does not work with single bonds, 
but with clusters of bonds. It has been noted earlier that 
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the distribution of load over the receptors in a finite-sized 
cluster induces non-trivial cooperativity between adhesion 
bonds: as one bond is disrupted, the load on the remain- 
ing bonds increases and cluster stability diminishes [2,9, 
[T0l[TTj. Likewise, as a new bond is established, the other 
bonds feel less force and cluster stability is enhanced. 

For many common receptor-ligand systems, the lig- 
ands are tethered to a surface through polymeric spac- 
ers. This reduces the disturbance of ligand structure to 
preserve its specificity and allows exploration of space 
for receptors so that the effective affinity of surfaces cov- 
ered with specific bonds is increased. For example, it has 
been shown recently in a macroscopic shearing experiment 
for streptavidin-biotin bonded beads that bonding is en- 
hanced if ligands are tethered with polymeric spacers [12] . 
Thus the distance between ligands and receptors is an im- 
portant determinant for specific adhesion. When cells ad- 
here to and eventually spread on a surface, the distance be- 
tween the ligand- and receptor-carrying surfaces decreases 
from the /im- to the nm-range on the time scale of minutes 
13] . The final receptor-ligand distance for integrins is of 
the order of 15 — 20 nm, which is bridged by the polymer 
spacer carrying the ligand. 

During recent years, the distance-dependent binding of 
tethered ligands has been investigated experimentally as 
well as theoretically. In experiments using the surface force 
apparatus and high-affinity streptavidin-biotin bonds [141 
115] . it has been shown that binding depends on rare, strongly 
extended conformations of the spacers and that the onset 
of binding is followed by a very fast increase of the fraction 
of bound tethers. Binding of tethered ligands to receptors 
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has been described using a combination of Monte Carlo 
simulations and reaction-diffusion theory. Simulations of 
pearl-bead chains confined between two walls were used 
to determine the force-extension relation for the polymers 
and to derive a potential landscape for the movement of 
the ligands. Later reversibel tethered bonds have been 
treated as deep but finite potential wells in the polymer 
potential landscape derived from simulations and used in 
the reaction-diffusion equations [161IT7] . Moreover, the ef- 
fect of changing receptor-ligand distance with a prescribed 
velocity has been discussed. It was concluded that kinetic 
effects are most important for strong binders with large 
affinity while thermodynamic equilibrium dominates for 
weak, reversible binders with small affinity. Recently, a 
theoretical treatment has also been given for receptor- 
ligand binding between curved surfaces, where different 
bonds are not equivalent for geometrical reasons [18j . 

It is important to note that these theoretical stud- 
ies have been conducted in the framework of a mean- 
field description for a large numbers of independent, non- 
cooperative bonds. No theoretical treatment has been pre- 
sented yet for the impact of receptor-ligand distance on 
the stability and dynamics of finite-sized clusters with co- 
operative bonds. Experimentally, it is well known that the 
lateral arrangement of the integrins and therefore the in- 
tegrin cluster size is regulated by cytoplasmic proteins like 
talin and a-actinin, which can bind both to the integrins 
and to the actin cytoskeleton [19] . Using image correlation 
microscopy, it has been shown for migrating cells that the 
integrins which are not yet organized in adhesions are al- 
ready preclustered with an average cluster size of three to 
four [20]. With a measured area density of few hundreds 
of integrins per [im 2 , this corresponds to a lateral distance 
well above 100 nm between the different mini-clusters. As 
the adhesion contacts nucleate and grow, integrins are in- 
creasingly clustered, until they approach a density of 1.000 

- 10.000 per /xm 2 , corresponding to a lateral distance of 10 

- 30 nm. The progression from very few to thousands of in- 
tegrins per cluster suggests that finite size effects might be 
highly relevant in the stabilization of initial cell adhesion. 

In this paper, we theoretically study the effect of feed- 
back and cooperativity on the distance-dependent rececepto: 
ligand binding dynamics in finite-sized adhesion clusters. 
To this end, we use a one-step master equation, which has 
been used before to study the adhesion cluster stability 
under mechanical load, both for constant force 1121 and 
linearly rising force [22j . Here we extend this framework 
to include the effect of receptor-ligand distance. We start 
in Sec. [2] by introducing the appropriate one-step master 
equation. In addition, we derive two corresponding contin- 
uum descriptions, namely a Fokker-Planck equation and a 
deterministic differential (mean field) equation. Our mod- 
eling framework can be used for any spatial distribution 
of the ligand in the direction normal to the substrate to 
which it is tethered. For simplicity, we start with a har- 
monic tether potential (spring model), which is the first 
order approximation for all polymer models at small ex- 
tensions. In Sec. [31 we analyze the stationary states of the 
mean field equation and derive one-parameter bifurcation 



diagrams which show that the receptor-ligand dynamics of 
adhesion clusters leads to bistability between bound and 
unbound states. In Sees. [5] and [6] we discuss stationary 
and dynamic properties of the master equation, respec- 
tively. In particular, we find that large adhesion clusters 
are stabilized in the bound state due to very large switch- 
ing times to the unbound state. In Sec. [7] we combine our 
conceptual framework with the worm-like chain model for 
semiflexible polymers to study the effect of finite poly- 
mer contour length. We conclude in Sec. [8] by discussing 
some biological applications of our results. A short report 
on our main results regarding the spring model has been 
given before [23] . 



2 Model 

2.1 Master equation and continuum limits 

Fig. [T] shows a schematic representation of our model for 
an adhesion cluster. The mechanical properties of the force 
transducer holding the receptors and of the polymeric 
tethers holding the ligands are represented by harmonic 
springs. We consider a situation in which N t receptor- 
ligand pairs are arranged in parallel between the planar 
force transducer and the substrate. At a given time, each 
of these adhesion bonds can either be open or closed. All 
bonds are considered to be equivalent so that the state 
of the adhesion cluster is characterized by the number of 
closed bonds i alone. The number of open bonds is Nt — i. 
Because i ranges from i = (completely unbound state) to 
i — N t (completely bound state) , our model has N t + 1 dis- 
crete states. Every bond changes its state (rebinds or rup- 
tures) through thermally activated, stochastic transitions. 
The stochastic rates for these transitions will be specified 
below. Therefore, the stochastic variable i changes by dis- 
crete steps ±1 (one-step process). The time dependence of 
the functions {piit)}^ representing the probability that 
i bonds are closed at time t are described by a one-step 
master equation which has the general form 

dpi 

-T^ = n+iPi+i + gi-iPi-i - {n + g l }p l . (1) 

The forward rates gi for the formation and the reverse 
rates for the rupture of a closed bond have to be speci- 
fied for all states < i < N t . They have the general form 
fi = ik ff(i) and gi = (Nt — i)k on (i). The reverse rate is 
the product of the number of closed bonds and the sin- 
gle bond off-rate k ff(i) because the closed bonds rupture 
independently and at a given time any of them could be 
the next to break. The forward rate is the product of the 
number of open bonds (Nt — i) and the single bond on- 
rate k on (i) because the open bonds are independent and 
any one of them could be the next to bind. Both, off- and 
on-rate will in general depend on i. The form of the tran- 
sition rates implies that the system is constrained to the 
interval < i < N t by two natural, reflecting boundaries 
at i — and i — N t , provided the single bond rates are 
finite at these boundaries. 



Thorsten Erdmann, Ulrich S. Schwarz: Impact of receptor-ligand distance on adhesion cluster stability 



3 




Fig. 1. (a) Cartoon of a cell adhering to a substrate through two sites of adhesion. Each adhesion corresponds to an elastic 
deformation in regard to the average cell-substrate distance, (b) Simple mechanical model for a site of adhesion. The cluster 
consists of Nt bonds (here N t = 5). At a given time t, i of these bonds are closed (here 3) while Nt — i (here 2) are open. 
The force transducer at the top (cell envelope) and the ligand tethers at the bottom are modeled as harmonic springs, with 
rest lengths It and lb and spring constants kt and kb, respectively. The equilibrium positions of the springs (dashed lines) are 
separated by the receptor-ligand distance I. The extensions of the bond springs and the transducer are denoted by x;, and 
x t = I — Xb, respectively. 



If the transition rates are continuous functions of i, 
that is Qi = g(i) and = r(i), it is useful to introduce 
a mean-field description for the average number of closed 
bonds, N(t) = (i) = 53i=i The time derivative of 

N(t) follows from the master equation Eq. fl} in terms of 
averages over the transition rates. Dragging the average 
into the argument, it can be approximated as 



dN 
~dt 



(S(i))-<r(i)) « ff«i»-r«»» = g(N)-r(N) . (2) 



This procedure yields an ordinary differential equation for 
N which in the following we term the mean field equation. 
It is an exact equation for the mean number of bonds N 
if the transition rates are linear functions in i. Otherwise, 
it is only valid in the limit of large system size Nt ■ 

To assess the role of fluctuations in the continuum 
limit, we derive a Fokker-Planck equation from the mas- 
ter equation. To achieve this, the one-step master equation 
Eq. ((T|) is written in operator form [24) 



dp{i,t) 
dt 



{(B-lMO + ^-l^jpfM) (3) 



where the operators E and E 1 act on the index i as 

Ef(i):=f(i+1) and E^f(i) := f(i - 1) . (4) 

In a continuum description for i these can be expanded in 
a Taylor series, 

oo oo 

E = l + J2-<d? and E^ 1 = 1 + ^ — r~9" . (5) 

n— 1 n=l 

Inserting Eq. ([5]) in Eq. ([3]) leads to the Kramers-Moyal 
expansion [25 



dt 



E 



,2n+l 



(2n+ 1)! 



A{i) 



oo o2n I 



in which the Kramers-Moyal coefficients are defined by 

A(i) = g(i) - r(i) and D(i) = g(i) + r(i) . (7) 

The Kramers-Moyal expansion Eq. (J6j) is an exact contin- 
uum representation of the master equation. Terminating 
the infinite series after second order yields the Fokker- 
Planck equation 

°EM = {-a i A(i) + ±%D®}p(i,t) (8) 

where A(i) is the drift and D(i) the diffusion coefficient, 
which determine the short time behavior of the first and 
second centered moments, respectively: 



and 



A(i(t )) = Jim 



D(i(t )) = jjm o 



(i(t + At) - i{t Q )) 
At 



({i(t + At) - i(t )) 2 
At 



(9) 



(10) 



(6) 



By deriving the differential equations for average and vari- 
ance directly from the master equation, one can verify that 
the left hand sides are the derivatives dN/dt and do 2 /dt 
for the initial conditions Nito) = i(to) and cr 2 (io) = at 
time t = to 24,21]. The Fokker-Planck equation thus de- 
scribes the average and the variance correctly. However, 
the centered moments of higher order are not correctly 
represented due to the truncation procedure. 

In general, a continuous description of a one-step pro- 
cess through the mean field equation Eq. @ or the Fokker- 
Planck equation Eq. © will be valid if the variation of 
transition rates and probability distribution is small over 
the step size Ai = ±1. Below we will use the mean field 
equation for a bifurcation analysis in order to demonstrate 
that our system is bistable. Because the Fokker-Planck 
equation results from a truncated Kramers-Moyal expan- 
sion, it cannot be used to describe the full dynamics. For a 
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bistable system it has been shown before that the diffusion 
approximation Eq. ([8]) leads to an error in the stationary 
probability distribution and an overestimation of the tran- 
sition rates between the coexisting states [26) . However, 
this error is small for sufficiently smooth transition rates. 
Moreoever, the extrema of the stationary distribution are 
correctly described by the diffusion approximation of the 
Fokker-Planck equation as used here. Therefore we will 
use it below to investigate how the stationary states are 
affected by thermal fluctuations which are not included in 
the mean field approximation. In order to describe the full 
stochastic dynamics, we will use the original master equa- 
tion Eq. ([1]). Alternatively, an improved Fokker-Planck de- 
scription could be used as explained in Ref. [26] . 



2.2 Transition rates 

2.2.1 Reverse rate 

The bond dissociation dynamics is characterized by the re- 
verse rate r(i) and strongly depends on the forces acting 
in the cluster. We start with the simplest possible model 
as suggested by the cartoon in Fig. []] that is the trans- 
ducer and ligand tethers are modeled as harmonic springs 
with rest lengths it and £b and spring constants k t and fcb, 
respectively. Because the receptors are usually firmly at- 
tached to the actin cytoskeleton, deformation of the trans- 
ducer requires a local deformation of the whole cell mem- 
brane as shown in the cartoon of Fig. QJa). Therefore the 
stiffness k t of the transducer is the combined stiffness of 
plasma membrane and cell cortex. The extensions Xb for 
a bound ligand and xt for the transducer satisfy the rela- 
tion Xb + x t = £, where I is the receptor-ligand distance in 
the completely dissociated state. In the following, we will 
treat the relaxed or unloaded receptor-ligand distance £ as 
a parameter which has been fixed externally, e.g. by the 
average cell-substrate distance or in an experiment with 
the surface forces apparatus. Mechanical equilibrium re- 
quires 



iF b = F t = k t x t = k t (£ - x b ) 



so that 



Fb(i) = hx b {i) 



k h £ 



Fb/h) (11) 



(12) 



l + i(k b /kt) 

The dissociation rate of a bond under force is described by 
the Bell model as fc Q // = k exp (Fb/F ), where ko is the 
unstressed off-rate and Fq is the bond's internal force scale 
[2] . The Bell model can be rationalized in the framework of 
Kramers' theory for the escape over a sharp energy barrier 
[3]. The reverse rate r(i) for the one-step master equation 
now follows as 



r(i) = ik exp (F b (i)/F ) 



(13) 



When a bond ruptures, force is redistributed according to 
Eq. (fT2|) over the smaller number of bonds so that the 
load on the remaining bonds increases. As a consequence, 
the extension of the remaining tethers increases which in- 
creases the force even further. A decrease of the number 



of closed bonds thus increases the off-rate for the single 
bonds and makes rupture of a further bond more likely. 
Thus the reverse rate from Eq. (fl"3|) describes a positive 
feedback mechanism for rupture which results from coop- 
erativity in load sharing. 

2.2.2 Forward rate 

The bond association dynamics is characterized by the for- 
ward rate g(i) and is strongly determined by the receptor- 
ligand distance which has to be bridged for the formation 
of a new bond. In the cartoon of Fig. [TJ the ligands are at- 
tached to springs and can move between the ligand-coated 
surface at the bottom and the transducer surface at the 
top. Hence, they move in a truncated harmonic potential, 
which is 

U(x) = ys 2 if -4 < x < x b (14) 

and U (x) — oo otherwise. For a finite temperature T, 
the probability density for the ligand to be at position 
x e [—tbi x b] reads 



p{x) = i exp {-U{x)/k B T) = -| exp (-k b x 2 /2k B T) 



with the partition sum 



(15) 



Z 



irksT 
2h 



erf 



_Kblb_ 

2k R T 



erf 



kbx\ 



2k R T 



(16) 

Here erf (a;) = (2/y / 7r) J Q exp (— t 2 ) dt is the error func- 
tion. The binding process can conceptually be divided into 
two steps. First, ligand and receptor have to come suffi- 
ciently close to form an encounter complex [27.28,29ll30|. 
Second, this entangled state has to react to form the fi- 
nal complex. For a stationary distribution of tethers, the 
first step is limited by p(xb), the probability of the lig- 
and to be close to the transducer surface. The second step 
is described by the on-rate k on for the case that ligand 
and receptor are sufficiently close within the binding ra- 
dius £ b ind- The forward rate g(i) for the one-step master 
equation Eq. (QJ thus reads 



g(i) = k on (N t - i)p{x b ) ■ 



(17) 



When an open bond rebinds, force is redistributed over 
the larger number of bonds according to Eq. (|12[) and the 
load on the single bonds decreases. This reduces the ex- 
tension of the bound tethers and increases the density of 
free ligands close to the transducer. An increase of the 
number of closed bonds thus increases the on-rate for the 
single bonds and makes rebinding of a further bond more 
likely. Thus the forward rate from Eq. (TTT)) describes a pos- 
itive feedback mechanism for binding which results from 
cooperativity in the formation of encounter complexes. 

2.3 Adiabatic assumption 

In the derivation of reverse and forward rates it was as- 
sumed that relaxation of transducer and polymer tethers 
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to mechanical equilibrium after a change in i is fast com- 
pared to rupture and rebinding of adhesion bonds. This 
assumption is a prerequisite for the definition of discrete 
states described by the number of closed bonds alone and 
hence for the validity of the master equation. Experimen- 
tally it has been found for a biomembrane force probe that 
the damping time for the system with tethered bonds is 
on the order of 10 -3 s [3T], which is typically at least 
one order of magnitude smaller than the transition times 
for adhesion bonds [32] • The polymer relaxation time is 
determined by the Zimm time [33] . For a Flory chain, it 
is tz — riRp/kBT where n is the viscosity of the sur- 
rounding fluid, Rf is the Flory radius Rf — aN - 6 with 
Kuhn length a and number of Kuhn segments N. The 
viscosity of water is on the order of 77 ~ 10~ 3 Pas. For 
polymers with Rp ~ 10 nm and for fc^T ~ 4 pNnm one 
has Tz ~ 10~ 7 s. This result agrees with earlier estimates 
for polyethylene glycol chains pT l UB l HH], The very fast 
relaxation time scale for the polymers allows us to use a 
stationary density distribution for the ligands. 



2.4 Dimensionless parameters 

The master equation Eq. {T} together with Eqs. (fT3|) and 
fT2"|) for the reverse rate and Eqs. fI7J). ([IS]) and (JTSJ) for 
the forward rate completely specifies our model. For the 
following, it is useful to introduce dimensionless quanti- 
ties. First we introduce dimensionless time r = k$t. Then 
we non-dimensionalizc all distances by writing them in 
units of the unstressed ligand tether length lb. The di- 
mensionless relaxed ligand-receptor distance is denoted 
by A = 1 1 lb- We also introduce k = kb/kt, the ratio of 
the spring constants of bonds and transducer, and <f> — 
kb£b/Fo, the force in units of Fq that is necessary to ex- 
tend a ligand tether spring by lb, i.e. to twice its unstressed 
length. Then, the reverse rate reads 



r(i) = i exp (4>\/ (1 + Ki)) — i exp (4>X(i)) 



(18) 



where X(i) = A/ (1+ni) has been introduced as an abbrevi- 
ation for the extension of the tethers (receptor-ligand dis 
tance with i closed bonds). Regarding the association pro- 
cess, we define the dimensionless on-rate 7 = (k on /ko)(£bi n d 
which is weighted with the ratio of binding radius lund 
and unstressed tether length lb, and the inverse thermal 
energy ksT non-dimensionalised by the tether energy at 
an extension equal to their rest length, (3 
Then, the forward rate reads 



k b e b /2k B T. 



g{i) = 2 7 (/V t ~ i) 



exp I 



-/3A 2 (i)) 



erf 



(l3i) +erf (pi\(i) 



(19) 



With the definition of these dimensionless rates, the dy- 
namical equations, that is master equation, mean field 
equation and Fokker-Planck equation have the same form 
as in Eqs. ([1}, ([2]) and (|SJ), but with time derivatives in r 
rather than t. Since reverse rate Eq. (|18p and forward rate 
Eq. (|19p are both non- linear in i, the mean field equation 
for N is only valid for large system size. 



Our model now contains six dimensionless parameters. 
The number of receptor-ligand pairs is given by the clus- 
ter size N t . The conditional rebinding rate 7 describes 
the rate of binding with a flat density distribution (infi- 
nite temperature) on an interval of length lb- The relative 
stiffness of the tethers, K — kb/kt, implies the two limits of 
k — > 00 (soft transducer) and n — > (stiff transducer). In 
the following we will use the intermediate case k = 1 . The 
dimensionless force constant <f> measures the force needed 
to stretch the tethers to twice their unstressed length in 
units of the intrinsic force scale of the adhesion bonds. 
For an entropic spring, this essentially scales as the ratio 
of two length scales, the bond reactive compliance ksT/Fo 
and the rest length of the tethers. In practice it will have a 
rather small value and in the following we use <f> = 0.1. The 
mobility of the ligands is represented by the inverse ligand 
temperature (3. In the limits of soft transducer (k — > 00) 
and very high ligand temperature (/3 — > 0), our model 
simplifies to a case which we have studied before in order 
to assess adhesion cluster stability under force [TTll21| . In 
this paper, we rather focus on the role of ligand-receptor 
distance A, which in combination with the different spring 
constants replaces the dimensionless force / used in the 
earlier model. In experimental setups, A is certainly the 
most accessible parameter. The parameters and their def- 
inition are summarized in Tab. I. There, we also give the 
typical range of parameters which was used for calcula- 
tions and some estimates for the integrin-fibronectin sys- 
tem. 



3 Bifurcation analysis of the mean field 
equation 

To analyze the stationary solutions of the mean field equa- 
tion, we first consider the dependence of the reverse rate 
r(N) from Eq. |TSJ and the forward rate g(N) from Eq. (fl"9"|) 
on the number of closed bonds N and on the model pa- 
rameters. Stationary solutions are the fixed points of the 
Eq. ([2]) and correspond to intersections of r(N) and g(N), 
because then the time derivative dN/dr = g(N) — r(N) 
/vanishes. Fig.^a) plots r(N) and g(N) as function of N 
for different values of the relaxed receptor-ligand distance 
A. The single bond off-rate, k off (N) = exp (0A/(1 + kN)), 
is finite at iV = 0, thus r(0) = 0. With increasing N, the 
reverse rate r(N) increases almost linearly, although the 
single bond off-rate is a monotonous decreasing function of 
N. The weak influence of the exponential off-rate k ff(N) 
is mainly due to the small force constant = 0.1: the en- 
tropic tether force does not suffice to accelerate bond rup- 
ture appreciably. For larger <f> the reverse rate increases 
quickly from r(0) = and can have a local maximum 
and minimum at small TV. For large N the linear term 
dominates in any case. Alternatively, non-monotonous be- 
havior could be induced by a large n, that is for a soft 
transducer. The forward rate g(N) vanishes for N — Nt, 
goes through a maximum at intermediate TV and then de- 
creases. At TV = the forward rates is always positive and 
approaches zero only in the limit of infinite f3 or A. At 
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definition 


typical 


integrins 


meaning 


P := k b l'i/2k B T 


0.1... 10 


7 


inverse ligand temperature in units of tether energy 


A := £/i b 


0.1... 10 


0.75 


receptor-ligand distance 


N t 


10... 25 


10 


cluster size 


<f> := k b l b /F 


0.1 


0.27 


tether force in units of internal force scale of bonds 


K :— kb/kt 


1 


0.9 


ratio of tether and transducer stiffness 


7 : = 7(4md/4) 


1 


1 


conditional single bond on-rate 



Table 1. The six parameters of the model: definitions, typical values used in the calculations, estimates for the integrin- 
fibronectin system and meaning. 




Fig. 2. (a) Reverse and forward rate r(N) and g(N) and (b) the time derivative dN/dr = g(N) — r(N) as a function of the 
number of closed bonds N for Nt = 10 and = 1. The receptor-ligand distance is A = 1, 2.5 and 5. The other parameters are 
7 = 1, k = 1 and (f> — 0.1. 



intermediate values for A, there are three intersections of 
g(N) and r(N). Fig. E^b) plots dN/dr as function of N 
for the same set of parameters as in Fig.[2la). It is positive 
at N = 0, because r(0) = and g(0) > 0. At intermediate 
N and small A, the time derivative has a maximum which 
reflects the maximum of g(N) before it becomes negative 
at large N where the reverse rate dominates. The fixed 
point at large N is stable and represents the bound state 
of the adhesion cluster containing a large number of closed 
bonds. For intermediate A there are three fixed points, in- 
cluding two stable fixed points at large N (bound state) 
and N w (unbound state). The two stable fixed points 
are separated by an unstable one. At large A there is only 
one stable, unbound state which approaches N = in the 
limit of large A. 

Fig. [3] summarizes the behavior of the fixed points in 
the form of two one-parameter bifurcation diagrams which 
show the fixed points as function of A and f3, respectively. 
At small separation, a single stable fixed point exists at 
large N. Here, the adhesion cluster is bound because the 
force on the bonds is small and the density of free ligands 
close to the receptor is large, therefore rupture events are 
rare and can be balanced by rebinding. With increasing 
A, the force on the bonds increases and the number of 
closed bonds in the bound state decreases. At large A, 
there is a single stable fixed point at JV w 0. Here, the 
adhesion cluster is unbound because forces are large and 
ligand density at the receptors is small so that rupture 
events occur frequently and cannot be balanced by re- 



binding. The transition from bound to unbound proceeds 
via two saddle-node bifurcations. At small A, the stable 
unbound fixed point appears together with an unstable 
fixed point separating the stable ones. The unstable fixed 
point merges with the bound stable fixed point at larger A. 
In the window of bistability between the two bifurcations, 
two stable fixed points coexist. The position of this win- 
dow moves to smaller A with increasing j3. The bifuraction 
behavior as function of the inverse ligand temperature (3 
is qualitatively similar although the position of the stable 
bound state initially increases with (3 because the forward 
rate initially increases. With decreasing A, the position of 
the bistable range shifts to larger (3 while its width in- 
creases strongly. 

In Fig. [4] we construct a stability diagram which iden- 
tifies the region of bistability as a function of A and (3. The 
positions of the bifurcation delineating the bistable region 
were determined numerically. For large (3 the width of the 
interval in A stays almost constant while the position de- 
creases slowly which explains the very large bistable range 
in (3 at small A. In general, A and (3 are inversely related. 
For the lower bifurcation one can use the approximate cri- 
terion that the slope of dN/dr with respect to N has to 
become negative. Neglecting r'(N), one is left with the 
condition dg(N)/dN = at N = 0. For large (3 > 1 one 
finds A 2 ~ 1/(2(3k). For (3 <C 1 and for sufficiently large A, 
on the other hand, one finds (3 ~ ln(KiV t )/A 2 . Thus in both 
limits, (3 and A are related by an inverse square root. At 
small (3, the width of the bistable range decreases as the 
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Fig. 4. Stability diagram for adhesion clusters: solid lines are numerically determined positions of the lower and upper bifurca- 
tions as function of /3 and A. The shaded area between the curves is the region of bistability, above this region there is a single, 
unbound state while below it a single stable bound cluster exists. The curves are derived for the cluster size N t = 10, k = 1, 
4> = 0.1 and (a) 7 = 1 and (b) 7 = 5. 



two curves eventually converge and bistability vanishes. 
The stability diagram can be regarded as the projection 
of the cusp-like surface of the fixed points on the (A,/3) 
plane. For large clusters and also for larger on-rates 7, the 
curves will meet at negative (3 and bistability persist for 
positive p. 

In the framework of biochemical control of biological 
systems, bistability is commonly associated with an un- 
derlying positive feedback mechanism [34] . In our case, 
bistability can arise from two positive feedback mecha- 
nisms as described above. First, there is positive feedback 
for bond rupture: as one bond breaks, the force on the 
remaining bonds increases, thus increasing their dissocia- 
tion rate. Second, there is positive feedback for binding: as 
one ligand binds a receptor, the receptor-ligand distance 
is decreased and the binding rate for the other ligands is 
increased. In general, we verified that in our model, both 
mechanisms can lead to bistability. However, for the pa- 
rameter range chosen here it is only the positive feedback 
of binding which is responsible for the observed bistability. 



As shown in Fig. [2j the reverse rate r(N) increases almost 
linearly for the set of parameters chosen. The forward rate 
g(N), on the other hand, is non-monotonous. Thus for the 
parameter range chosen here, the positive feedback under- 
lying bistability is mostly due to the forward rate g(N). 

4 Bifurcation analysis with a stochastic 
potential 

4.1 Stationary solution of the Fokker-Planck equation 

The Fokker-Planck equation Eq. (|8]) has the stationary 
solution 

p -v=wM 2 iw> di ')- <2o) 

where C is a normalization constant. The integrand in the 
exponent exists and the expression is integrable because 
the Fokker-Planck coefficients are finite and defined on a 
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compact interval. In the absence of sources and sinks the 
flux 

J*(i) = A(t)P s (i) - \ {D'{i)P s {i) + D(i)P*{i)>} (21) 
has to vanish. The derivative of P s (i) with respect to i is 

p^r- m L D ' {t) p s (^ (22) 



D(i) 



so that indeed 



J*(i) = - \ {D'(i) + 2A(i) - D'(i)}\ P s (t) = . 

(23) 

In particular, the flux through the boundaries vanishes, as 
required for reflecting boundaries. 



4.2 Stochastic potential 

The stationary probability distribution Eq. (|20|) can be 
used to define an energy landscape E(i) by 



E(i) 



log P s (i). 



(24) 



The extrema of this potential are determined by the con- 
dition 



dE(i) 
di 



1 dP a 



P s (i) di 



= & 



dP s (i) 
di 



= 0. (25) 



With Eq. (f22|) for the first derivative of the distribution, 
the position of the extrema of the stochastic potential E(i) 
are thus determined by 



A(i)--D'(i) 



0. 



(26) 



These extrema have a similar physical meaning as the 
fixed points of the mean field equation, but they are more 
rigorous in including thermal fluctuations. In the frame- 
work of the stochastic potential, bistability requires a bi- 
modal potential landscape in which two minima of the 
stochastic potential coexist. The coexisting minima (max- 
ima of the probability distribution) are separated by a 
potential barrier (minimum of the probability distribu- 
tion) . The separated peaks of the probability distribution 
in these minima are commonly referred to as macrostates 
of the stochastic system because they are possible real- 
izations of a macroscopic, deterministic system. One can 
regard the extrema as the fixed points of the dynamical 
system 

j t = A(i) - \D'(i) = g(i) - r(i) - \{g'{i) + r'{i)) . (27) 

This is the mean field equation corrected for the effects 
of non-homogeneous mobility. Eq. (|2"T|) allows to deter- 
mine extrema in the same way as the fixed points of the 
mean field equation in the previous section. For a con- 
stant diffusion coefficient the fixed points are identical to 
the extrema of the stochastic potential. Non-homogeneous 
diffusion terms change the position of fixed point and can 
even destroy fixed points or create new ones [noise-induced 
transitions) [35] . 



4.3 Bifurcation analysis 

Using the stochastic dynamic system, the topology of the 
stochastic potential can be analyzed just as the mean field 
equation by deriving bifurcation diagrams showing the po- 
sitions of the macrostates (the extrema of the potential) 
as function of the model parameters. The extrema are de- 
termined as stable and unstable fixed points of Eq. (j2"7| . 
Fig.[5]shows these stochastic bifurcation diagrams as func- 
tion of receptor-ligand distance A in comparison to those 
from the mean field equation. For sufficiently small A, 
the upper stable fixed point of the deterministic equation 
agrees well with a bound macrostate at large N in the 
stochastic potential. The saddle-node bifurcation in which 
the bound macrostate vanishes occurs at much smaller 
separation A than in the deterministic picture, thus fluc- 
tuations destabilize the adhesion cluster. The minimum of 
the stochastic potential lies above the unstable determin- 
istic fixed point and becomes negative at small A. For pos- 
itive i, the stochastic potential has no second minimum, 
but it has a boundary minimum at i = 0. This second 
macrostate corresponds to the unbound fixed point of the 
mean field equation. For larger clusters, the agreement be- 
tween the fixed points of the stochastic potential and the 
mean field equation improves. In the limit of very large 
clusters, the two solution approach each other as shown 
in Fig. [S] for N t = 100 and 250. In this case, the unbound 
state is practically indistinguishable from i = 0. 



5 Stationary solutions of the master equation 

For a one-step master equation on a finite range without 
sources or sinks, stationarity Pi(oo) = implies detailed 
balance, that is r(i)pi(co) = g(i— l)pi_i(oo). Iterating this 
relation results in the stationary probability distribution 

^-4^1144 ^0<i<N t . (28) 
po(oo) r{i) £X r{j) 

The normalization constant is the stationary state proba- 
bility for the completely dissociated state i = 0, 



I , , f(0) TT 3{j) 



(29) 



Fig. [7] shows a density plot of the stationary distribution 
{pi(oo)}^* as function of relaxed receptor-ligand distance 
A for cluster sizes N t — 10 and N t = 25. For small A, 
there is a single peak at a finite number of closed bonds 
which is broadened by fluctuations. This corresponds to 
the bound state of adhesion clusters. For large A, there 
is a single maximum at the completely dissociated state 
i = 0, which is the unbound state. In an intermediate 
range of A, the stationary distribution has two maxima 
and bound and unbound adhesion clusters coexist. Thus 
the full stochastic model indeed shows bistability as sug- 
gested by the bifurcation analysis of the mean field equa- 
tions. For the smaller cluster size Nt = 10, fluctuations are 
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Fig. 5. Comparison of the extrema of the stochastic potential with the deterministic fixed points as function of (a) A with 
P — 0.5, 1 and 2 and (b) f3 with A = 1, 1.5 and 2 for the cluster size Nt = 10. For small parameters A and /3, stochastic and 
deterministic results agree well for the stable fixpoints at large N, but differ strongly in the transition regions. 




Fig. 6. Extrema of the stochastic potential as function of A compared with the deterministic fixed points for the same parameters 
as in Fig. [5] but for cluster sizes (a) Nt = 100 and (b) Nt = 250. For these large system sizes, stochastic and deterministic 
results agree well over the full range of parameters. 



large and the transition from bound to unbound appears 
rather smooth. For larger systems, the transition becomes 
sharper and discontinuous. This discontinuity is demon- 
strated by the average number of closed bonds in the adhe- 
sion cluster which shows a steep decrease as the occupancy 
switches from bound to unbound. For a bimodal distribu- 
tion with two distinct macrostates, average numbers of 
closed bonds can be defined in the two peaks separately. 
We use the probability functions fe(oo)}^* 3 for the upper 
and {pj(oo)}*_ for the lower peak with proper normal- 
ization to calculate the average number of closed bonds 
in the two macrostates. Fig. [7] shows that these averages 
vary slightly with A and the steep decrease of the full aver- 
age is mostly due to the change in occupancy probability 
than in the position of the peaks. Fig. [7] also shows the 
bifurcation result from the mean field equation. For the 
larger system N t = 25, the position of the maxima are 
in good agreement with the fixed points and the onset of 
bistability at small A agrees well with the first bifurcation. 
For the smaller system Nt = 10, the agreement between 
fixed points and maxima is still good, but the onset of bi- 



modality is overestimated by the lower bifurcation. Here, 
the stochastic potential yields a much better estimate for 
the locations of the bifurcation (not shown). For growing 
cluster size, the range of parameters in which the coexist- 
ing macrostates are occupied to a similar degree shrinks, 
although the range of bistability as revealed by the de- 
terministic equation or the stochastic potential grows. In 
analogy to first order phase transitions, the discontinuous 
transition from bound to unbound will occur at a sharp 
parameter value in the limit of infinite clusters. 

Biological systems are usually of finite size and time 
dependent processes are often important. Whether coex- 
isting states and transitions between them can be observed 
depends on the relative occupancy of the coexisting states 
and on the time scale for transitions between the states. 
The bare presence of bistability identified in the determin- 
istic equation is irrelevant if the time scale for the tran- 
sitions is larger than the typical observation times. The 
time scale for transitions will also be important if param- 
eters are time dependent. If the change of parameters is 
faster than the time scale for equilibration of the probabil- 
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Fig. 7. Density plot of the stationary probability distribution of the master equation {pi(oo)}°^ as function of separation A 
for cluster sizes (a) N t = 10 and (b) Nt = 25. The other parameters are /3 = 1.0, 7 = 1, k = 1 and (f> = 0.1. Dark regions 
indicate high probability. The curves are the average number of closed bonds in the bound and unbound macrostates, in the 
full distribution and as predicted from the mean field equation, respectively. 



ity distribution, mctastable states will be populated and 
hysteresis in the transition parameters will be observed. In 
the following section, we therefore analyze dynamic prop- 
erties of the master equation. 



6 Dynamic properties of the master equation 

The stationary probability distribution Eq. arises by 
averaging over many individual trajectories or over a sin- 
gle trajectory for a long time. Fig. [8] shows two sample 
trajectories for N t = 10 bonds and = 0.5 at A = 1 and 
A = 3, respectively. The stochastic trajectories are gener- 
ated with the Gillespie algorithm for exact stochastic sim- 
ulations 36,37] and show how the number of closed bonds 
i changes over time. For A = 1, the stationary distribu- 
tion is unimodal and has a single peak around the bound 
state. Due to the small cluster size, the number of closed 
bonds fluctuates strongly around the average; occasionally 
it reaches the completely dissociated state, but rebinding 
takes place immediately. This leads to a single, broad peak 
as in Fig. [7] at small A. For A = 3 the stationary distribu- 
tion is bimodal. The sample trajectory alternates between 
bound and unbound state. As long as the time spent in the 
respective states is large enough, the time taken for the 
actual transition is negligible. When bound, the trajectory 
fluctuates around an average as in (a). An encounter of the 
completely dissociated state, however, is usually followed 
by a longer time with no or few bonds. Single closed bonds 
are formed occasionally, but this rarely leads to the for- 
mation of a large number of closed bonds. Increasing the 
receptor-ligand distance further increases the time spent 
in the unbound state relative to that in the bound state. 
This means that the occupancy probability of the bound 
state is reduced. The trajectories alternating between the 
states with very short transition times yield the bimodal 
stationary distribution. The probability to find the system 
in one of the states is proportional to the time spent in 



that state before a transition. For very large A with uni- 
modal distribution, binding of trajectories does not take 
place with appreciable frequency. 

The dynamic properties of the stochastic model can 
be characterized by the mean first passage time T m ^ n be- 
tween two states m and n, that is, the time it takes on 
average to reach n for the first time from m. To elucidate 
the relevance of the fixed points in the stochastic system, 
the transition times from the unbound state m = to the 
bound state n ~ N t /2 and vice versa are of particular in- 
terest. The value N t /2 can be used because in this context 
it is only relevant that the bound state is above the tran- 
sition barrier between the two macrostates; the dynamics 
within the respective bassins of attraction are much faster 
than the dynamics across the barrier. For a one-step mas- 
ter equation like Eq. {1} , the mean first passage time from 
the initial state m to the final state n satisfies the recur- 
sion relation [53] 

g(n) {Tm+i,™ - r„,m} + r(m) {T m _i,„ - T m ,„} = -1 

.(30) 

with the boundary condition T m ^ m — 0. For a transition 
from a state to to n > to, that is for an increase of the 
number of closed bonds, one has 

T ~-t{m+'im&M- (31) 

For the reverse transition from a state to to n < m, where 
the number of closed bonds decreases, one has 

The first term in curly brackets in Eqs. (|3"Tj) and (|32p is 
the mean first passage time for a trajectory exclusively 
with binding or rupture, respectively. The second term 
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Fig. 9. Mean first passage times T m ,n from Eqs. (|3ip and (|32[) for transitions between unbound state m — and unbound state 
n — N t /2 as function of A for (a) Nt = 10 and /3 = 0.5, 1 and 2 and (b) for f3 = 1 with N t = 10, 25 and 50. The other parameters 
are 7 = 1, k = 1 and (f> — 0.1. 



with the product over the ratio of rebinding and rupture 
rates describes the increase of the mean first passage time 
through backward reactions, that is rupture if to < n and 
rebinding if to > n. 

Fig. [9] plots the mean first passage times T m , n from 
Eqs. (j3Tj) and (|32|) for transitions from the unbound state 
to = to the bound state n = N t /2 (binding time) 
and from the bound state to the unbound (unbinding 
time). For parameter values for which only a single sta- 
ble macrostate exists, the transition time into this state 
is very small and on the order of magnitude of transi- 



tions between neighboring states. The time for transitions 
in the reverse direction becomes extremely large for large 
clusters. The range of bistability is characterized by bind- 
ing and unbinding times which are both larger than the 
single step transition times. Because barrier crossing itself 
is a fast process, the ratio of binding to unbinding time 
equals the ratio of occupancy of the two macrostates in 
the range of bistability. The point where the transition 
times are equal thus defines a stochastic transition point 
at which both states are equally occupied. The plots of the 
transition times as function of A for different f3 at N t — 10 
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in Fig. [9^a) show that this transition point shifts to larger 
A with decreasing /3, because with increasing temperature 
(ligand mobility), larger separations can be bridged by 
the ligand tethers. Fig. [9^b) plots the transition times for 
j3 = 1 at different N t . This shows that the stochastic tran- 
sition point shifts to larger separation values with growing 
cluster size. At the same time the corresponding switch- 
ing times grow super-exponentially fast. This implies that 
for a given separation and growing cluster size, the bound 
macrostate will effectively become the only stable state; 
frequent switching between bound and unbound states can 
thus only occur for small clusters. It is important to note 
that these important conclusions can only be drawn by 
considering the full stochastic dynamics. 

Stochasticity of binding and strong dependence of the 
mean first passage time on A has to be considered in mea- 
surements of the binding ranges of tethered ligands. If the 
receptor-ligand distance is reduced step by step and the 
transducer is held at a constant distance during a short 
period of time t s , binding will typically be observed at a 
distance where the mean first passage time from unbound 
to bound T(X) is comparable to r s . In general, binding 
is always possible for all distances smaller than the con- 
tour length of the tethers. The actual binding distance is 
a stochastic variable. If transitions proceed with the con- 
stant rate 1/T(A), the probability to observe binding after 
the n th step at a distance A„ is 

n-l 

p = (1 - cxp (-r s /T(A„))) [J exp (-T s /T(Xi)) . (33) 

i=l 

With the strong dependence of the binding time on A as 
shown in Fig. [9] this distribution will be have a sharp peak 
when T(X n ) < r s < T(Aj<„). The probability to bind 
at any distance below a given A„ is then close to a step 
function as observed experimentally and theoretically [15] 

cog. 

7 Extension to worm-like chain model 

Unlike in the simple harmonic spring model, real polymers 
are not infinitely extensible, but are characterized by a 
finite contour length. A commonly used model for real 
polymers is the so-called worm-like chain or Kratky-Porod 
model [53]. It has been used before to model semi- flexible 
biopolymers like DNA, F-actin or titin [38 . A worm- like 
chain is characterized by the contour length L and the 
persistence length L p which describes the bending stiffness 
of the filament. The worm- like chain model can also be 
extended to include elasticity of the polymer backbone 
which allows stretching beyond the contour length [39] . 
The forces needed to stretch the polymer monomers are 
much larger than the typical thermal forces and will not 
be considered in the following. 

7.1 Force extension relation and rupture rate 

For a worm-like chain, the force F w i c which induces an 
average extension x of the worm-like chain polymer can 



be approximated by the interpolation formula [38] 

„ . . knT f 1 x I 1 , 

With the first term in curly brackets the force diverges as 
x approaches the contour length L. The second term de- 
scribes harmonic behavior at small extensions with a force 
constant 3ksT '/2LL p . The third, constant term guaran- 
tees that the force vanishes at vanishing extension. We 
express Eq. (|34p in non-dimensional units by writing ex- 
tension in units of the contour length, £ = x/L, and force 
in units of the intrinsic force scale Fq of adhesion bonds, 
fwic — F w i c /Fq. The force extension relation Eq. ([34]) then 
reads 

fwUO = 4>wic 1 4( - 1 ]_ ^ 2 + £ - ^ j (35) 

where the ratio 4> w i c — fcsT/(FoL p ) is defined in analogy 
to </>. For small extension, 1, the constant of propor- 
tionality between £ and f w t c is 3c/>/2. 

The extension of bound tethers is determined by 
mechanical equilibrium between tethers and transducer. 
In non-dimensional units this reads 

^W=W + « 4 >-i}-*---«'> < 36 > 

which has to be solved for The parameter X w i c — i/L 
is the non-dimensional relaxed receptor-ligand distance 
and k w i c — (3kBT/2LL p )/kt measures the ratio of the 
harmonic force constant of the polymer and the force con- 
stant of the transducer. The two parameters are analogous 
to A and k for the spring model. Solving Eq. ([36]) for 
yields tether extension and force f w lc(£b(i)) as function of 
the number of bound tethers alone. This result has to be 
inserted into Bell's expression fc Q // = e^ wlc ^', leading to 
the reverse rate of the adhesion cluster 

r(i) = ie /ro!cW (37) 

which has to be used in the one-step master equation 
Eq. ([1]) and the mean field equation Eq. J2}. As a poly- 
nomial of third order, Eq. ([36]) can be easily solved for 

&(*)■ 

7.2 Rebinding rate 

The energy needed to stretch a worm-like chain from zero 
extension to an extension Xb is given by the integral over 
the force F w i c 

V wlc (x) = f ^ } F wlc {x')dx' . (38) 
Jo 

In non-dimensional units the energy is 

« (em- 1 [ 1 I 6(0 n 
v wlc (^))- Tp { A(1 _ m) + - 2 r - l ) 

(39) 
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where the dimensionless persistence length is A p = L p /L. 
With the Boltzmann factor e~ Vwlc ^ the density of un- 
bound ligands at the transducer is 

pfa(i)) = e-*»«"«*W»/Z(&(i)) (40) 

where 

Z(Sb(i)) = / e^df + / e -««i.tt)df (41) 
Jo Jo 

is used as the partition sum. The first term is added to 
prevent the density from diverging at £h = 0; for a 
harmonic potential with the same force constant as for 
the worm-like chain was used. It takes into account that 
for entropic reasons the ligands are found on average at a 
certain height above the substrate. The exact distribution 
of the polymers in the presence of a wall is unknown [40] . 
In fTS] , it has been calculated by Monte Carlo simulations 
for a bead-pearl model. It was found that at large exten- 
sions, the distribution resembled that of a freely jointed 
chain [33] . The exact form of the distribution should have 
no influence on the generic aspects we are interested in. 
The distribution below the rest length is used only for 
normalization but is irrelevant for binding. 

The forward rate of the adhesion cluster as function of 
the number of closed bonds i is 

g(i) = j(N t - i)p(i) , (42) 

where 7 = {k on /k )(£bind/L) has been defined. Eqs. ([37]) 
and (|4"2"]) together with the master equation Eq. ([T]) and 
the deterministic equation Eq. ([2]) define the dynamical 
system. The six dimensionless parameters for the worm- 
like chain model are defined in analogy to those for the 
harmonic model. In the following we use A p = 1, 7 = 1, 
Kwic = 1-5 and 4> w i c = 0.1. The choices for k w i c and <fi w i c 
mean that the properties of the tethers at small extension 
are similar as above for the harmonic tethers. 



7.3 Analysis of the worm-like chain 

A steady state analysis of the mean field equation Eq. ([2J) 
with the transition rates Eq. ([37]) and Eq. l[4"2"]) from the 
worm-like chain model shows that the bifurcation behaviour 
is very similar to the one obtained for the linear springs, 
except that distances larger than X w i c = 1 are not possible 
due to the finite contour length. Fig. [10] shows a density 
plot of the stationary probability distribution Eq. ([28]) for 
the worm-like chain transition rates Eq. ([37]) and Eq. (|42[1 
as function of X w i c . Together with the distribution, the de- 
pendence of the upper stable fixed point of the mean field 
equation and the average number of closed bonds of the 
full distribution as well as that of the bound and unbound 
macro-states are displayed. Again the binding region is 
bounded by the maximum extension X w i c — IjL = 1. 
As for the harmonic tethers a bimodal region is found in 
which two macrostates coexist. The average number of 
closed bonds in the bound state and the position of the 



maximum agree well with the upper stable fixed point 
of the deterministic equation. Both depend hardly on the 
receptor-ligand distance X w i c . The average number of closed 
bonds jumps from bound to unbound state in a discontinu- 
ous transition. This transition becomes sharper for larger 
clusters, that is, the width of the region in which both 
states are occupied to an appreciable degree decreases 
with increasing cluster size. 

The physical reason for the striking plateau in the 
bound state is the non-linearity of the worm-like chain 
force extension relation, which reflects the strain-stiffening 
typical for biopolymers. A binding tether thus pulls the 
relatively soft transducer until the stiffnesses match, that 
is to the regime of harmonic tethers with n w ic = 1.5. 
Therefore, the final extension of the bound tethers de- 
creases less than linear with X w i c and the effect on the 
number of closed bonds is weaker than for the harmonic 
model. If the transducer was replaced by another worm- 
like chain, the behavior of the system should be more like 
the harmonic model. 



8 Conclusion 

In this paper, we have introduced a stochastic model which 
allows to study the interplay of cooperative binding and 
unbinding for finite-sized adhesion clusters mediated by 
tethered ligands. Our model is based on established prin- 
ciples of receptor-ligand binding, including Kramers-type 
rupture rates and separation-dependent binding rates based 
on the notion of an encounter complex. By implement- 
ing these principles in the framework of a one-step mas- 
ter equation, we were able to apply many powerful tech- 
niques from stochastic dynamics, including a Kramers- 
Moyal derivation of a Fokker-Planck equation (which in 
turn corresponds to a stochastic potential) and exact solu- 
tions for stationary solutions and mean first passage times. 
In particular, a bifurcation analysis based on the stochas- 
tic potential could be compared to the bifurcation analysis 
of the mean field equation to the master equation. 

Our model shows that the simple mechanical model of 
Fig. [T] can lead to a bistable situation in which two differ- 
ent states of adhesion, bound and unbound, coexist. The 
underlying reason for the occurence of bistability is the 
existence of two mechanisms for positive feedback, one for 
rupture and one for binding. Cooperative bonds share the 
force exerted by the transducer so that the force exerted 
on each closed bond reduces upon binding. At the same 
time, the extension of the tethers reduces upon binding, 
which then increases the probability for further binding. 
In consequence, the transition rates in the master equa- 
tion are both strongly non-linear functions of the number 
of closed bonds, which both can lead to positive feedback. 
This leads to an instability for the intermediate numbers of 
closed bonds and thus to bistability. The model discussed 
in this paper is an extension to a previously introduced 
model in which only the rupture rate was a non-linear 
function of the number of closed bonds [TT1|21| . It could 
be re-obtained from the current model in the limit of a 
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Fig. 10. Stationary probability distribution {pi(oo)}^ from Eq. (|28|) with the transition rates Eq. (|37[) and Eq. (|42|l for the 
worm-like chain model and for cluster sizes (a) Nt = 10 and (b) Nt = 25 plotted as a function of the relaxed receptor-ligand 
distance X w i c . The other parameters are k w i c = 1.5, <f)wic = 0.1, A p = 1 and 7=1. The curves show the upper stable fixed point 
of the mean field equation, the average number of closed bonds in the bound and the unbound macrostate, and in the whole 
cluster, respectively. 



soft transducer, that is for k — ► oo. In the present dis- 
cussion, non-linear effects were mainly due to cooperative 
rebinding while the non-linearity of the rupture rate was 
weak. 

In the mean field description, bistability leads to hys- 
teresis for the binding and unbinding range: binding from 
the unbound state takes place at a smaller distance than 
unbinding from the bound state. Unlike for previous dis- 
cussions for non-cooperative bonds [16|[T7] , this is not due 
to kinetic effects. In the stochastic description, the bistable 
system is characterized by a bimodal stationary probabil- 
ity distribution. Fluctuations over the barrier separating 
the stable states allow equilibration of the distribution. 
The time scale for equilibration has been calculated as 
the mean first passage time between the different adhe- 
sion states. Due to these finite transition times kinetic 
effects from changes of external parameters are important 
for the behavior of the system. In analogy to thermody- 
namic systems at first order phase transitions, changes 
of parameters on time scales smaller than the transition 
times allow occupation of the metastable states, while for 
slower changes, a stationary distribution will occur. In the 
bimodal region, the system frequently alternates between 
bound and unbound state. For large systems, the transi- 
tion times become very large and the transition between 
bound and unbound becomes very sharp. The bimodal re- 
gion in which two states are occupied to an appreciable 
and comparable amount becomes very small. 

The predictions of this paper for the internal dynam- 
ics of adhesion clusters can be investigated with exper- 
imental setups like AFM or the surface force apparatus 
which allow to study the specific binding of two opposing 
surfaces with controled separation. However, to map out 
the occupancy distribution for bound and unbound states 
could be tedious because of large transition times and low 
occupancy outside the dominant state. If measuring the 
binding range by stepwise reduction of the receptor-ligand 



distance, binding will take place when the transition time 
is smaller or equal to the waiting time. For large systems 
this will lead to hysteresis as in the deterministic case. For 
smaller systems, hysteresis will decrease due to the smaller 
transition times. In experiments with the surface forces 
apparatus [T4l[T5] , a behavior similar to the one predicted 
here has been observed. After initial binding a large at- 
tractive force was measured until an equilibrium position 
of the surfaces was established. In our model this increase 
would be due to the increased extension of the transducer 
spring. In those experiments, irreversible bonds with very 
large affinity have been used so that repeated transitions 
between bound and unbound states could not be observed. 

In order to test our results, it would be necessary to 
use reversible adhesion bonds with low affinity. In general, 
binding through reversible bonds is highly relevant for bi- 
ological systems. In particular, cell-matrix adhesion is me- 
diated by reversible bonds like the ones between integrin- 
receptors in the plasma membrane and fibronectin in the 
extracellular environment. In this case, the equilibrium 
length of the ligand tether is lb = 11 ran. Using the pa- 
rameter values given in Tab. I, one then finds from the 
stochastic model for N t = 5 that for small adhesion clus- 
ters, bistability should occur around a relaxed receptor- 
ligand distance of 8 nm (A = 0.75) [23]. Together with 
the ligand and receptor rest lengths, this results in a cell 
substrate distance of around 20 nm, that is the physio- 
logical value for cell-matrix adhesion [13]. Therefore the 
mechanism of bistability as described here can be used by 
cells to explore the extracellular space by many small and 
transient adhesions. On encountering favorable conditions, 
these small adhesions might mature, e.g. by recruitment 
of additional receptors. The results presented here show 
that this quickly leads to switching times which keep the 
adhesions in the bound macrostate. 

In order to present a simple and reasonable model, 
here we have made the crucial assumption that all bonds 
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are equivalent. This assumption leads to a one-step master 
equation and allows the application of many powerful tools 
from stochastic dynamics. In the future, our model could 
be extended to include additional aspects of biomimetic 
or biological systems, which usually however can not be 
described in the framework of a one-step master equa- 
tion. In order to describe the growth and shrinkage of 
adhesion clusters, the overall number of bonds N t should 
be made dynamic, possible involving regulation through 
the cytoskeleton. Assuming an elastic rather than a rigid 
transducer requires solution of elastic equations in order 
to derive the exact details of the force distribution. The 
force distribution would also be changed when accounting 
for possible curvature of the opposing surfaces. Finally it 
would also be interesting to consider the effect of disor- 
der, e.g. in bond resting lengths or single bond on- and 
off-rates. 
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